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Interpocket polarization model for magnetic structures in rare-earth hexaborides 
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The origin of peculiar magnetic structures in cubic rare-earth (R) hexaborides RBg is traced 
back to their characteristic band structure. The three sphere-hke Fermi surfaces induce inter- 
pocket polarization of the conduction band as a part of a RKKY-type interaction. It is shown 
for the free-electron-like model that the interpocket polarization gives rise to a broad maximum 
in the intersite interaction I{q) around q = (1/4, 1/4, 1/2) in the Brillouin zone. This maxi- 
mum is consistent with the superstructure observed in R=Ce, Gd and Dy. The wave-number 
dependence of I{q) is independently extracted from analysis of the spin-wave spectrum mea- 
sured for NdBe. It is found that I{q) obtained from fitting the data has a similarly to that 
derived by the interpocket polarization model, except that the absolute maximum now occurs 
at (0,0, 1/2) in consistency with the A-type structure. The overall shape of I{q) gives a hint 
toward understanding an incommensurate structure in PrBg as well. 
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Interest in cubic rare-earth (R) hexaborides RBg 
comes mainly from the rich structure in their ordered 
phases. The best studied example is CeBg which under- 
goes a quadrupole (orbital) order at 3.3 K and then a 
magnetic order at 2.4 K.^^ The magnetic ground state 
is characterized by double-fc structure with wave vectors 
(1/4, ±1/4, 1/2) in units of the reciprocal lattice param- 
eter 27r/a. Since the orbital order is superimposed on 
the magnetic order, it has been suspected that the or- 
bital degeneracy in the crystalline electric field (CEF) 
ground state Fg plays an important role.^-* Recently, 
however, neutron scattering experiment on GdBg has de- 
tected an equivalent wave vector in the ordered phase 
below 15 K.'^) Since the trivalent Gd ion has a half-filled 
4f shell without orbital degrees of freedom, the order at 
k = (1/4,1/4,1/2) should have an origin which does 
not depend so much on the particular configuration of 4f 
electrons. It is known that PrBg also has the same wave 
number in the magnetically ordered ground state below 
4.2 K, but the intermediate phase between 4.2 K and 
6.9 K has an incommensurate structure.*^ On the other 
hand, the ground state of NdBg has a simple antifer- 
romagnetic structure called the type I (or A-type) with 
alternating plane polarized along and against (0,0,1).^'^^ 

In this paper we propose a simple model to under- 
stand the origin of these structures from a unified point 
of view. The basic observation is that the Fermi surface 
of RBg consists of three nearly spherical pieces centered 
on the X points = (1/2, 0, 0), = (0,1/2,0),X2 = 
(0, 0, 1/2) in the Brillouin zone. The RKKY interaction 
involves interpocket polarization, which has a new char- 
acteristic wave vector = (1/2, 1/2, 0) which connects 
Xx and Xy, and equivalent ones. Just like the ordinary 
RKKY interaction can bring about the antiferromag- 
netic ordering by halving the reciprocal lattice vector, 
the halving of the characteristic wave vector (1/2, 1/2,0) 
can bring about the ordering at (1/4, 1/4, 1/2). 

Let us consider the case of GdBg where the 4f electrons 
have only the spin degrees of freedom. The exchange 
interaction between a 4f-spin Si at Ri and conduction 



electrons is taken to be 
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where N is the number of lattice sites, and JWk,^ 
determined by the exchange integral involving 4f wave 
functions and the conduction states. We follow the pre- 
vious argument^-' to derive JWk,p- In analogy with the 
APW method we consider a muffin-tin sphere centered at 
the origin. The Bloch function ifikif) of the conduction 
band is expanded inside the sphere as 
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where Rk\ (r) describes the radial part with orbital index 
A, and Y-^!^\r) with f — r/r is the cubic harmonics for 
the point-group representation F and its component 7. 
We neglect the /c-dependence of Rk\{r) since the extent 
of 4f electrons is smaller than that of 5d electrons which 
contribute dominantly to the exchange. Because the or- 
bital angular momentum is zero in Gd'^"'", the exchange 
integral becomes diagonal with respect to the azimuthal 
quantum number of 4f states, and to (F,7). Thus the 
exchange interaction becomes isotropic with a factor 
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The RKKY interaction I{q) is given by 

, I •■i'fc,fe-|-ql 
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with /(e) being the Fermi function. The intrapocket con- 
tribution to I{q) comes from such k and k + q that be- 
long to the same pocket of the conduction band. In addi- 
tion, there arises the interpocket contribution explained 
earlier. Let us take the free-electron-like dispersion and 
set |Wfc,fc_|_qp constant in order to see the consequence 
of the interpocket contribution in the simplest manner. 
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In order to keep the lattice periodicity, we take sum- 
mation over the reciprocal lattice vectors G rather than 
restricting the fc-summation within the Brillouin zone. 
The G-summation corresponds to inclusion of higher en- 
ergy bands. Namely we introduce 

ii{q)=Y,F{q + G)n{q + G), 
G 

where F{q + G) is a form factor to be specified later, and 
n(q) is the Lindhard function multiplied by the partial 
density of states at the Fermi level. The Fermi wave num- 
ber kp is given by kpa/ir = n~^/^ = 0.9656/\/2, which 
means that the three spherical Fermi surfaces barely 
touch with one another. In the real RBg system, the 
Fermi surface also contains a fine structure along (1,1,0) 
and equivalent directions.^' 

Adding both intrapocket and interpocket contribu- 
tions we obtain for q G Brillouin zone: 

3 

I{q) = J2 [3U{q) + 2 ^ II(q - K,)] = j\{q), 

i=l 

where 311 (q) accounts for the three equivalent pock- 
ets, and II(q — Ki) describes the interpocket polariza- 
tion. The K,'s are given by Ki = (0, 1/2, 1/2), = 
(1/2, 0,1/2), Kg = (1/2,1/2,0) in units of 27r/a. The 
factor 2 for the interpocket term enters because each 
pocket can be both starting and ending states of the tran- 
sition. For simplicity we take the form factor such that 
F{k) = 1 if \ka\ < Qir/a for all components a = x,y,z 
and zero otherwise. The choice of the cut-off in the form 
factor hardly influences the q-dependence of the RKKY 
interaction, although it does influence the absolute value. 
Specifically with F(q + G) — 1 for all G, II(q) diverges 
logarithmically by summation over G. 

Figure 1 shows x{q) in the X-M-R plane of the Bril- 
louin zone with X = (0,0, 1/2), M = (1/2,0,1/2) and 
R = (1/2,1/2,1/2). The unit of ordinate is such that 
n(0) ~ 1, and the large numerical value of x(q) comes 
from summation over G. For the intcrsite interaction, 
only the variation in the q-space is relevant since the 
average of x{q) represents the intra-site contribution. 

We have also made a scan of x{q) afong (1/4, 1/4, q^) 
and found that the peak indeed occurs at = 1/2. It 
is apparent that the interaction favors the magnetic or- 
der near the center of the X-M-R plane, namely around 
(1/4,1/4,1/2). Since the ridge extends more toward 
M rather than i?, an incommensurate structure with 
qx ^ qy can be realized by slight change of the system 
parameters. 

We now analyze in more detail the character of the 
conduction band, which consists mainly of t2u molecular 
orbitals of 2p electrons in Bg clusters hybridized with Cg 
orbitals of 5d electrons. One of the t2u orbitals has the 
angular dependence z(a;^ — y^) if seen from the center of 
the Bg cluster, and hybridizes best with the 5d — 
orbital at neighboring rare-earth sites. Since z{x'^ — y"^) 
changes sign below and above the B4 plane, the wave 
number (0,0, 1/2) gains the bonding energy optimally.^^ 
Thus the bottom of the conduction band goes to X^ and 
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Fig. 1. Wave-number dependence of the intersite interaction in the 
X-M-R plane: (a) three-dimensional plot; (b) contour plot. 



equivalent points. For F = Cg we use a simplified no- 
tation (7|fc) = (i[?^''(fe) with 7 being either the state 
x"^ — y^ or — . Then we have a large amplitude 
(a;^ — y'^\Xz), while {iz^ — r^\Xz) is negligible. At an- 
other point Xx — (1/2,0,0), the wave function has the 
character of y^ — which can also be represented by 

\y'-z')^-l\x'-y')-^\3z'-r'), (4) 

with use of the basis set at Xz- Thus the orbital flip 
from \x'^ — y"^) to |t/^ — z^) can take place even with the 
cubic symmetry. We note that the finite overlap does 
not contradict with the orthogonality of Bloch functions 
with different k. 

The relative weight of the interpocket polarization 
against the intrapocket one should influence the detailed 
behavior of /(q). We estimate from the above argument 
the weight factor Wk,p for k = Xz and p = Xx relative 
to Wk,k as 

Wk,plWk,k - -1/2. 

These points Xz and Xx, however, are not on the Fermi 
surface. In the region where two pieces of the Fermi sur- 
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face almost touch with each other, the iuterpockct con- 
tribution connecting the nearby k states should have a 
larger weight factor than the intrapocket one with remote 
k and p. We have made a tight-binding calculation tak- 
ing the eg and t2u orbitals, and examined the character 
of the wave functions at various points in the Brillouin 
zone. It is found that 2p-electron weight is larger than 
the 5d-electron weight in general, and the latter changes 
gradually from \x'^ — y"^) to jy^ — z^) as fe moves from 
Xz to Xx- In a future work, we shall evaluate Wk,p by 
using realistic wave functions. 

The presence of orbital degeneracy in rare-earth ions 
other than Gd makes it necessary to consider more com- 
plicated form of H^f. Namely not only the spin exchange 
interaction but multipole interactions also enter.'^'^^'^^^ 
As long as the conduction band consists purely of Cg for 
the 5d electron part, the wave-number dependence of the 
multipole intersite interactions is the same as that of the 
exchange interaction. Actually, however, t2g{= Fs) also 
enters into eq.(2). In the presence of orbital degener- 
acy, Wk,p is no longer given by eq.(3) but with different 
weights for each Moreover, hybridization between 4f 
electrons and boron 2p electrons may become important 
in the open-shell case. The hybridization constitutes an- 
other mechanism of the intersite interaction. 

With these complications in mind we proceed to anal- 
ysis of the exchange interaction in NdBe where the spin- 
wave spcctnim has been measured. The dipole part of 
H^f can be taken in the same form as eq.(l) except that 
Si is replaced by the angular momentum operator J, 
with J = 9/2. With only the magnetic intersite inter- 
action, the easy axis of the magnetic moment should be 
along (1, 1, 1).^^^ Actually the moment is parallel or an- 
tiparallel to (0,0, 1), which has been explained in terms 
of ferroquadrupolar interaction.^^' With inclusion of 
the magnetic and Fs-type quadrupolar interactions, we 
consider the following model: 

where we assume that the average of the quadrupole mo- 
ment does not depend on a site. Other interactions such 
as the Fs-type quadrupolar interaction^^^ are neglected 
since they do not affect the spin-wave spectrum. 

The CEF ground state is Fg^-*, which is four-fold de- 
generate, and the first excited state lies 132~135 K 
above. "'^^^ We introduce the Pauli matrix a" (a = 
X, y, z) to describe the Kramers pair, and another Pauli 

(2) 

matrix r" to describe the orbital pair in the Fg quar- 
tet. Then the angular momentum operator J" within 
the Fg ' subspace is written as 

where ^ = —0.883 and rj = —4.712 are numerical con- 
stants corresponding to the Lea-Leask-Wolf^^^ parame- 
ter X = —0.82.^'^) The orbital effect on the magnetic 
moment is described by T" with = and 

j'X,y ^ _^zi2 ± V3t^/2. 



In the Neel state the Fg quartet undergoes a Zee- 
man splitting by the molecular field. This splitting in- 
duces a finite quadrupole moment which is enhanced by 
positive ^g. Then the lowest level is characterized by 
(t^,(7^) = (-|-,T) in the A-sublayer and (+, i) in the B- 
sublayer. One may expect two branches corresponding 
to excitations 

(+,ct) ^ (±,-ct), 

with intensities /±. Here cr =t,i in the A- and B- 
sublayers, respectively. The intensity ratio is 
given by 

/+//_ = (2C/ry- 1)73 ^0.13. 

Thus we identify the observed branch as the inter-orbital 
transition {+,cr) — > (— ,— a), and assume that intra- 
orbital branch was not detected because of the small 
intensity. 

By neglecting the small matrix elements for intra- 
orbital magnetic excitations, we obtain a reduced model 
which keeps only the two levels leading to /_. Assum- 
ing the A-type antiferromagnetic structure with Q = 
(0, 0, 1), we obtain the excitation spectrum by the stan- 
dard spin- wave theory as 

a;2 = [J(g)-A][J(q + Q)-A], (5) 

where J{q) is the Fourier transform of = (377^/16) 7^, 
and 

^=^-^^^J{Q) + l{i-n-^f9's. (6) 

The spectrum given by eq.(5) was also postulated in the 
previous work, where LOq experimentally measiircd was 
fitted by intersite interactions up to third neighbors. The 
authors of ref.6 noted that the calculated Neel tempera- 
ture was about half of the experimental one, 8^9 K. We 
point out further that the substantial softening around 
(1/4, 1/4,0) was not reproduced by the previous fit. 

In view of the fact that the RKKY interaction has a 
long range, we have included a sufficient number (= 34) 
of intersite interactions Jij. The comparison between 
theory and experiment is shown in Fig. 2. The experi- 
mental spectrum is well reproduced by our fit. The fit 
gives J{Q) = 1.16 meV and A = 1.47 meV. Then we ob- 
tain (/g = 108 mK from cq.(6). This value is in excellent 
agreement with experimental one, g'^ ~100 mK, deduced 
from elastic constant. In the mean-field approxima- 
tion the Neel temperature is not influenced by g'^, and is 
given by 

which is 9.3 K with the fitted value of J{Q)- 

Figure 3 shows J{q) in the X-M-R plane. The maxi- 
mum of J{q) occurs at (0, 0, 1/2) in consistency with the 
A-type order. In addition, there appears a local max- 
imum near (1/4,1/4,1/2). The latter indicates a ten- 
dency toward the ordering with q = (1/4, 1/4, 1/2), and 
brings about the softening of LOq near this wave num- 
ber. We note that the overall behavior of J{q), and thus 
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Fig. 2. The magnetic excitation spectrum in NdBg along various 
symmetry directions. Experimental data^' are shown by solid 
circles with error bars, and the solid lines are theoretical fits. 




M 



Fig. 3. Wave-number dependence of J{q) in the X-M-R plane. 



I{q) = (16/3?7)^ J(q), is similar to the result of the inter- 
pocket polarization model shown in Fig.l. This similar- 
ity supports relevance of the model to real RBe systems. 
The difference should mainly come from our simplifica- 
tion for Wk,p and e^, and partly from the presence of 
orbital degeneracy, hybridization, and correlation effect 
among conduction electrons. 

In this paper we have concentrated on the q- 
dependence of the intersite interaction I{q). As the 
magnetization grows, the associated nonlinearity favors 
a commensurate structure in general. Then the max- 
imum of I{q) does not necessarily give the ordering 
wave number at the ground state. We suggest that the 
incommensurate-commensurate transition in PrBg may 
be interpreted along this line. It should be worth inves- 
tigating detailed features which depend on 4f-clcctron 
configurations of each rare-earth species. 

Another feature to be addressed with finite order pa- 
rameters is the direction of magnetic moment at each 
site. Even with the same q = (1/4,1/4,1/2), the mo- 
ment patterns are rather different between CeBg^^ and 
DyBg.^^' While in CeBg the nearest-neighbor moments 
are orthogonal to each other and within the (001) plane, 
the moments in DyBg point to (1/2, 1/2, 1/2). The latter 
is consistent with the magnetic anisotropy in the para- 
magnetic region. In this connection it is interesting to 
inquire into the spin patterns of GdBg at low tempera- 
ture, since the magnetic anisotropy in the paramagnetic 
region is extremely small. 
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